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Abstract. We model magnetars as hybrid stars, which have a core of quark matter surrounded 
by hadronic matter. For this purpose, we use an extended version of the SU(3) non-linear 
realization of the sigma model in which the degrees of freedom change naturally from hadrons 
to quarks as the temperature/density increases. The presence of a variable magnetic field allows 
us to study in detail the influence of Landau quantization and the anomalous magnetic moment 
on the particle population of the star, more precisely on particles with different spin projections. 
This allows us to calculate the polarization of the system throughout different phases of the star, 
hadronic, quark and also a mixed phase. 



1. Introduction and Model Description 

Magnetars are compact stars that have extremely high magnetic fields. The highest magnetic 
field observed on the surface of a star is on the order of 10 15 G. The highest possible magnetic 
field in the center of stars, on the other hand, can only be estimated through models, even when 
applying the virial theorem. Some results indicate limiting magnetic fields ranging between 
B = 10 18 — 10 20 G [1-8]. Here we assume that magnetars are not necessarily composed 
of hadronic matter and describe them using a model that contains both hadronic and quark 
degrees of freedom. 

Within this approach, the hadrons (whole baryon octet) are replaced by quarks (up, down, 
strange) at high densities and/or temperatures. This happens as the effective masses of the 
hadrons increase and the effective masses of the quarks decrease within these limits. The 
aforementioned model (see Ref. [9] for details) is an extended version of the SU(3) non-linear 
realization of the sigma model within the mean- field approximation. Changes in the order 
parameters of the model a and $ signal chiral symmetry restoration and quark deconfinement, 
respectively. The potential for $ is an extension of the Polyakov loop potential [10] modified 
to also depend on baryon chemical potential. In this way our model is able to describe the 
entire QCD phase diagram, even at zero temperature. Fig. 1 shows that the model is in 
good agreement with lattice QCD constraints [11] and that it reproduces the liquid-gas phase 




Figure 1. (Color online) QCD Phase 
Diagram - Temperature as a function of 
baryon chemical potential showing first-order 
phase transitions. The dots represent critical 
points. 



Figure 2. (Color online) Particle densities 
function of baryon chemical potential as- 
suming global charge neutrality and chemical 
equilibrium at T=0. Quark densities are di- 
vided by 3. 



transition for symmetric matter. In this figure we also show results for neutron star matter, 
which is charge neutral and in chemical equilibrium. The phase transitions at low temperatures 
and high densities are of first order, whereas at high temperatures and low densities the model 
exhibits smooth crossovers. 

The SU(3) non-linear realization of the sigma model and its extension (that also contains 
quarks) have been successful in reproducing nuclear matter properties [12], heavy ion collision 
data [13], compact star and proto-neutron star properties [14-16]. As compact stars have 
temperatures of the order of 1 MeV, we can safely set their temperature to zero. As already 
mentioned, for star calculations we have to take into account charge neutrality and chemical 
equilibrium. Here we assume that the surface tension between the hadronic and quark phases 
is small [17] and allow charge neutrality to be global (only the combination of both phases sum 
up to zero charge). As a consequence, a mixed phase appears in the star. This can be seen in 
Fig. 2, which also shows that hyperons are almost completely suppressed by the appearance of 
the quarks. 

We include in the model a magnetic field in the z-direction that has varying magnitude. This 
is a more realistic approach than considering a constant magnetic field throughout the star and 
can prevent the creation of hydrodynamical instabilities due to pressure anisotropy [18-22]. This 
happens because, in our approach, the magnetic field only becomes extremely high in the center 
of the star, where the matter pressure is also high (see Ref. [23] for more details). More precisely, 
we assume an effective magnetic field B* that increases with chemical potential, running from 
a surface value B sur f = 69.25 MeV 2 = 10 15 G (when [Ib = 938 MeV) to different central values 
B c at large values of baryon chemical potential following [16] 

B*(n B ) = B surf + B c [l-e b 938 ], (1) 

with a = 2.5, b = —4.08 x 10 -4 and \ib given in MeV. As can be seen in Fig. 3, the values of 
the effective magnetic field only approach B c at very high baryon chemical potentials and, in 
practice, only about 70% of B c can be reached inside stars. The use of an explicit dependence of 
B on the baryon chemical potential instead of on density was chosen to prevent discontinuities 
in the magnetic field at the phase transition, where the baryon density is discontinuous. The 
constants a and b and the form of the B* expression were chosen to reproduce (in the absence 
of quarks) the effective magnetic field curve as a function of density from Rcfs. [5, 24, 25]. 
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Figure 3. (Color online) Effective magnetic 
field as a function of baryon chemical 
potential shown for different central magnetic 
fields. Recall that using Gaussian natural 
units 1 MeV 2 = 1.44 x 10 13 G. 



Figure 4. (Color online) Particle densities 
function of baryon chemical potential 



for B r 



5 x 10 5 MeV 



7.22 x 10 18 G 



including AMM. Black, blue and red stand 
for negative, while green, purple and orange 
stand for positive spin projections. Quark 
densities are divided by 3. 



2. Results and Conclusions 

The magnetic field in the z-direction forces the eigenstates in the x and y directions of charged 
particles to be quantized into Landau levels. The energy levels of all baryons are further split 
due to the alignment /anti-alignment of their spins with the magnetic field (anomalous magnetic 
moment effect, AMM). But even when the AMM is not taken into account, like in the quark phase 
in our model, only one of the spin projections contributes to the zeroth Landau level, creating 
a spin projection asymmetry in the system. In this work, we focus on the analysis of magnetic 
field effects on the chemical composition of the neutron star, the total spin polarization and the 
magnetization of the system. Studies of magnetic field effects on compact star observables can 
be found in Refs. [16,26-28]. 

The particle population is shown in Fig. 4 when a central magnetic field B c = 5 x 10 5 
MeV 2 = 7.22 x 10 18 G with AMM is considered. The "wiggles" in the charged particle densities 
mark the [is values, for which their Fermi energies cross the discrete threshold of a Landau level. 
The charged particle population is enhanced due to B, as their chemical potentials increase. 
Although the AMM is known to make the EOS stiffer, it does not have a very significant effect 
in the particle population. This fact can be easily understood in terms of polarization, when, 
instead of looking at the total particle density (sum of spin up and down particle densities) for 
each species, we look at the spin up/spin down particle densities separately. In this case some 
of these particles are enhanced while others are suppressed (see Ref. [29] for details). 

The total polarization of the system is defined as 

J2iQBi(p+i ~ P-i) 

F,iQBi(p+i + p-iY 

where Qb% stands for the baryon number of each species (3 for baryons, 1 for quarks and for 
leptons). We can see in Fig. 5 that the total polarization of the system (taking into account 
hadrons and quarks) is larger for larger magnetic fields and increases with chemical potential 
until the quarks appear. After this point the polarization decreases. This happens because 
the AMM for the quarks is not included in our calculations, as it is not fully understood for 
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Figure 5. (Color online)Total spin polariza- 
tion as a function of baryon chemical poten- 
tial for different central magnetic fields (with 
1 MeV 2 = 1.44 x 10 13 G) including AMM. 



Figure 6. (Color online) Effective magnetic 
field times magnetization as a function 
of baryon chemical potential for different 
central magnetic fields (with 1 MeV 2 = 
1.44 x 10 13 G) including AMM. 



these particles. In this case, the spin asymmetry comes exclusively from the spin asymmetric 
contribution to the zeroth Landau level. For an example of a complete study of the polarization 
without the AMM see Ref. [30]. The "wiggles" in Fig. 5 mark again when the Fermi energies of 
the charged particles cross the discrete threshold of a Landau level. 

Finally, we turn our attention to the calculation of the magnetization of the system. It is 
important to notice that in the case with the AMM, not only the magnetization of the charged 
particles has to be included, but also the magnetization of the uncharged particles. For details 
on the calculation of the magnetization including the AMM see Ref. [31]. Fig. 6 shows that the 
magnetization is larger for larger magnetic fields and it increases with baryon chemical potential. 
As in Figs. 4 and 5, the "wiggles" mark when the Fermi energies of the charged particles cross 
the discrete threshold of a Landau level. The magnetization of the system is a very important 
quantity as it relates with the pressure anisotropy of the system [18-22]. 

We have shown in this work some possible effects of strong magnetic fields in hybrid stars. 
The presence of different hadronic and quark degrees of freedom makes this quark-hadron sigma 
model an ideal tool for such an analysis in the different possible phases of the star. More 
specifically, we analyzed the effects of strong chemical potential dependent magnetic fields on 
particles with different spin projections, their total polarization and magnetization. 
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